WikiPathways 2024: next generation pathway database.

WikiPathways (wikipathways.org) is an open-source biological pathway database. Collaboration and open science are pivotal to the success of WikiPathways. Here we highlight the continuing efforts supporting WikiPathways, content growth and collaboration among pathway researchers. As an evolving database, there is a growing need for WikiPathways to address and overcome technical challenges. In this direction, WikiPathways has undergone major restructuring, enabling a renewed approach for sharing and curating pathway knowledge, thus providing stability for the future of community pathway curation. The website has been redesigned to improve and enhance user experience. This next generation of WikiPathways continues to support existing features while improving maintainability of the database and facilitating community input by providing new functionality and leveraging automation.

Human Monocytes Exposed to SARS-CoV-2 Display Features of Innate Immune Memory Producing High Levels of CXCL10 upon Restimulation.

INTRODUCTION: A role for innate immune memory in protection during COVID-19 infection or vaccination has been recently reported. However, no study so far has shown whether the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) can train innate immune cells. The aim of this study was to investigate whether this virus can induce trained immunity in human monocytes. METHODS: Monocytes were exposed to inactivated SARS-CoV-2 (iSARS-CoV-2) for 24 h, followed by a resting period in the medium only and a secondary stimulation on day 6 after which the cytokine/chemokine and transcriptomic profiles were determined. RESULTS: Compared to untrained cells, the iSARS-CoV-2-trained monocytes secreted significantly higher levels of IL-6, TNF-α, CXCL10, CXCL9, and CXCL11 upon restimulation. Transcriptome analysis of iSARS-CoV-2-trained monocytes revealed increased expression of several inflammatory genes. As epigenetic and metabolic modifications are hallmarks of trained immunity, we analyzed the expression of genes related to these processes. Findings indicate that indeed SARS-CoV-2-trained monocytes show changes in the expression of genes involved in metabolic pathways including the tricarboxylic acid cycle, amino acid metabolism, and the expression of several epigenetic regulator genes. Using epigenetic inhibitors that block histone methyl and acetyltransferases, we observed that the capacity of monocytes to be trained by iSARS-CoV-2 was abolished. CONCLUSION: Overall, our findings indicate that iSARS-CoV-2 can induce properties associated with trained immunity in human monocytes. These results contribute to the knowledge required for improving vaccination strategies to prevent infectious diseases.

Exploring pathway interactions to detect molecular mechanisms of disease: 22q11.2 deletion syndrome.

BACKGROUND: 22q11.2 Deletion Syndrome (22q11DS) is a genetic disorder characterized by the deletion of adjacent genes at a location specified as q11.2 of chromosome 22, resulting in an array of clinical phenotypes including autistic spectrum disorder, schizophrenia, congenital heart defects, and immune deficiency. Many characteristics of the disorder are known, such as the phenotypic variability of the disease and the biological processes associated with it; however, the exact and systemic molecular mechanisms between the deleted area and its resulting clinical phenotypic expression, for example that of neuropsychiatric diseases, are not yet fully understood. RESULTS: Using previously published transcriptomics data (GEO:GSE59216), we constructed two datasets: one set compares 22q11DS patients experiencing neuropsychiatric diseases versus healthy controls, and the other set 22q11DS patients without neuropsychiatric diseases versus healthy controls. We modified and applied the pathway interaction method, originally proposed by Kelder et al. (2011), on a network created using the WikiPathways pathway repository and the STRING protein-protein interaction database. We identified genes and biological processes that were exclusively associated with the development of neuropsychiatric diseases among the 22q11DS patients. Compared with the 22q11DS patients without neuropsychiatric diseases, patients experiencing neuropsychiatric diseases showed significant overrepresentation of regulated genes involving the natural killer cell function and the PI3K/Akt signalling pathway, with affected genes being closely associated with downregulation of CRK like proto-oncogene adaptor protein. Both the pathway interaction and the pathway overrepresentation analysis observed the disruption of the same biological processes, even though the exact lists of genes collected by the two methods were different. CONCLUSIONS: Using the pathway interaction method, we were able to detect a molecular network that could possibly explain the development of neuropsychiatric diseases among the 22q11DS patients. This way, our method was able to complement the pathway overrepresentation analysis, by filling the knowledge gaps on how the affected pathways are linked to the original deletion on chromosome 22. We expect our pathway interaction method could be used for problems with similar contexts, where complex genetic mechanisms need to be identified to explain the resulting phenotypic plasticity.

Elucidating the Corneal Endothelial Cell Proliferation Capacity through an Interspecies Transcriptome Comparison.

The regenerative capacity of corneal endothelial cells (CECs) differs between species; in bigger mammals, CECs are arrested in a non-proliferative state. Damage to these cells can compromise their function causing corneal opacity. Corneal transplantation is the current treatment for the recovery of clear eyesight, but the donor tissue demand is higher than the availability and there is a need to develop novel treatments. Interestingly, rabbit CECs retain a high proliferative profile and can repopulate the endothelium. There is a lack of fundamental knowledge to explain these differences. Gaining information on their transcriptomic variances could allow the identification of CEC proliferation drivers. In this study, human, sheep, and rabbit CECs are analyzed at the transcriptomic level. To understand the differences across each species, a pipeline for the analysis of pathways with different activities is generated. The results reveal that 52 pathways have different activity when comparing species with non-proliferative CECs (human and sheep) to species with proliferative CECs (rabbit). The results show that Notch and TGF-ß pathways have increased activity in species with non-proliferative CECs, which might be associated with their low proliferation. Overall, this study illustrates transcriptomic pathway-level differences that can provide leads to develop novel therapies to regenerate the corneal endothelium.

Allele-specific expression analysis for complex genetic phenotypes applied to a unique dilated cardiomyopathy cohort.

Allele-specific expression (ASE) analysis detects the relative abundance of alleles at heterozygous loci as a proxy for cis-regulatory variation, which affects the personal transcriptome and proteome. This study describes the development and application of an ASE analysis pipeline on a unique cohort of 87 well phenotyped and RNA sequenced patients from the Maastricht Cardiomyopathy Registry with dilated cardiomyopathy (DCM), a complex genetic disorder with a remaining gap in explained heritability. Regulatory processes for which ASE is a proxy might explain this gap. We found an overrepresentation of known DCM-associated genes among the significant results across the cohort. In addition, we were able to find genes of interest that have not been associated with DCM through conventional methods such as genome-wide association or differential gene expression studies. The pipeline offers RNA sequencing data processing, individual and population level ASE analyses as well as group comparisons and several intuitive visualizations such as Manhattan plots and protein-protein interaction networks. With this pipeline, we found evidence supporting the case that cis-regulatory variation contributes to the phenotypic heterogeneity of DCM. Additionally, our results highlight that ASE analysis offers an additional layer to conventional genomic and transcriptomic analyses for candidate gene identification and biological insight.

Guiding the choice of informatics software and tools for lipidomics research applications.

Progress in mass spectrometry lipidomics has led to a rapid proliferation of studies across biology and biomedicine. These generate extremely large raw datasets requiring sophisticated solutions to support automated data processing. To address this, numerous software tools have been developed and tailored for specific tasks. However, for researchers, deciding which approach best suits their application relies on ad hoc testing, which is inefficient and time consuming. Here we first review the data processing pipeline, summarizing the scope of available tools. Next, to support researchers, LIPID MAPS provides an interactive online portal listing open-access tools with a graphical user interface. This guides users towards appropriate solutions within major areas in data processing, including (1) lipid-oriented databases, (2) mass spectrometry data repositories, (3) analysis of targeted lipidomics datasets, (4) lipid identification and (5) quantification from untargeted lipidomics datasets, (6) statistical analysis and visualization, and (7) data integration solutions. Detailed descriptions of functions and requirements are provided to guide customized data analysis workflows.

Drug-target identification in COVID-19 disease mechanisms using computational systems biology approaches.

Introduction: The COVID-19 Disease Map project is a large-scale community effort uniting 277 scientists from 130 Institutions around the globe. We use high-quality, mechanistic content describing SARS-CoV-2-host interactions and develop interoperable bioinformatic pipelines for novel target identification and drug repurposing. Methods: Extensive community work allowed an impressive step forward in building interfaces between Systems Biology tools and platforms. Our framework can link biomolecules from omics data analysis and computational modelling to dysregulated pathways in a cell-, tissue- or patient-specific manner. Drug repurposing using text mining and AI-assisted analysis identified potential drugs, chemicals and microRNAs that could target the identified key factors. Results: Results revealed drugs already tested for anti-COVID-19 efficacy, providing a mechanistic context for their mode of action, and drugs already in clinical trials for treating other diseases, never tested against COVID-19. Discussion: The key advance is that the proposed framework is versatile and expandable, offering a significant upgrade in the arsenal for virus-host interactions and other complex pathologies.

Tissue-specific pathway activities: A retrospective analysis in COVID-19 patients.

The ACE2 receptors essential for SARS-CoV-2 infections are expressed not only in the lung but also in many other tissues in the human body. To better understand the disease mechanisms and progression, it is essential to understand how the virus affects and alters molecular pathways in the different affected tissues. In this study, we mapped the proteomics data obtained from Nie X. et al. (2021) to the pathway models of the COVID-19 Disease Map project and WikiPathways. The differences in pathway activities between COVID-19 and non-COVID-19 patients were calculated using the Wilcoxon test. As a result, 46% (5,235) of the detected proteins were found to be present in at least one pathway. Only a few pathways were altered in multiple tissues. As an example, the Kinin-Kallikrein pathway, an important inflammation regulatory pathway, was found to be less active in the lung, spleen, testis, and thyroid. We can confirm previously reported changes in COVID-19 patients such as the change in cholesterol, linolenic acid, and arachidonic acid metabolism, complement, and coagulation pathways in most tissues. Of all the tissues, we found the thyroid to be the organ with the most changed pathways. In this tissue, lipid pathways, energy pathways, and many COVID-19 specific pathways such as RAS and bradykinin pathways, thrombosis, and anticoagulation have altered activities in COVID-19 patients. Concluding, our results highlight the systemic nature of COVID-19 and the effect on other tissues besides the lung.

Understanding signaling and metabolic paths using semantified and harmonized information about biological interactions.

To grasp the complexity of biological processes, the biological knowledge is often translated into schematic diagrams of, for example, signalling and metabolic pathways. These pathway diagrams describe relevant connections between biological entities and incorporate domain knowledge in a visual format making it easier for humans to interpret. Still, these diagrams can be represented in machine readable formats, as done in the KEGG, Reactome, and WikiPathways databases. However, while humans are good at interpreting the message of the creators of diagrams, algorithms struggle when the diversity in drawing approaches increases. WikiPathways supports multiple drawing styles which need harmonizing to offer semantically enriched access. Particularly challenging, here, are the interactions between the biological entities that underlie the biological causality. These interactions provide information about the biological process (metabolic conversion, inhibition, etc.), the direction, and the participating entities. Availability of the interactions in a semantic and harmonized format is essential for searching the full network of biological interactions. We here study how the graphically-modelled biological knowledge in diagrams can be semantified and harmonized, and exemplify how the resulting data is used to programmatically answer biological questions. We find that we can translate graphically modelled knowledge to a sufficient degree into a semantic model and discuss some of the current limitations. We then use this to show that reproducible notebooks can be used to explore up- and downstream targets of MECP2 and to analyse the sphingolipid metabolism. Our results demonstrate that most of the graphical biological knowledge from WikiPathways is modelled into the semantic layer with the semantic information intact and connectivity information preserved. Being able to evaluate how biological elements affect each other is useful and allows, for example, the identification of up or downstream targets that will have a similar effect when modified.

A systematically derived overview of the non-ubiquitous pathways and genes that define the molecular and genetic signature of the healthy trabecular meshwork.

PURPOSE: The trabecular meshwork (TM) is situated in the most frontal part of the eye and is thought to play an important role in the regulation of the eye pressure. However, this tissue is rather difficult to harvest for research. The purpose of this study is therefore to integrate the existing gene expression data of the healthy TM to increase sample size and identify its signature genes and pathways. This provides a robust reference for the study of molecular disease processes and supports the selection of candidate target genes for new treatments. METHODS: A systematic search identified microarray data of healthy TM tissue. After quality control, datasets of low quality and deviating samples were excluded. Remaining individuals were jointly normalized and integrated into one database. The average gene expression of each tested gene over all individuals was calculated. The 25% genes with the highest average expression were identified as the most active genes in the healthy TM and used as input for pathway and network analysis. Additionally, ubiquitous pathways and genes were identified and excluded from the results. Lastly, we identified genes which are likely to be TM-specific. RESULTS: The gene expression data of 44 individuals, obtained from 18 datasets, were jointly normalized. Ubiquitous genes (n = 688) and ubiquitous pathways (n = 73) were identified and excluded. Following, 1882 genes and 211 pathways were identified as the signature genes and pathways of the healthy TM. Pathway analysis revealed multiple molecular processes of which some were already known to be active in the TM, for example extracellular matrix and elastic fiber formation. Forty-six candidate TM-specific genes were identified. These consist mainly of pseudogenes or novel transcripts of which the function is unknown. CONCLUSIONS: In this comprehensive meta-analysis we identified non-ubiquitous genes and pathways that form the signature of the functioning of the healthy TM. Additionally, 46 candidate TM-specific genes were identified. This method can also be used for other tissues that are difficult to obtain for study.

Systems Biology in ELIXIR: modelling in the spotlight.

In this white paper, we describe the founding of a new ELIXIR Community - the Systems Biology Community - and its proposed future contributions to both ELIXIR and the broader community of systems biologists in Europe and worldwide. The Community believes that the infrastructure aspects of systems biology - databases, (modelling) tools and standards development, as well as training and access to cloud infrastructure - are not only appropriate components of the ELIXIR infrastructure, but will prove key components of ELIXIR's future support of advanced biological applications and personalised medicine. By way of a series of meetings, the Community identified seven key areas for its future activities, reflecting both future needs and previous and current activities within ELIXIR Platforms and Communities. These are: overcoming barriers to the wider uptake of systems biology; linking new and existing data to systems biology models; interoperability of systems biology resources; further development and embedding of systems medicine; provisioning of modelling as a service; building and coordinating capacity building and training resources; and supporting industrial embedding of systems biology. A set of objectives for the Community has been identified under four main headline areas: Standardisation and Interoperability, Technology, Capacity Building and Training, and Industrial Embedding. These are grouped into short-term (3-year), mid-term (6-year) and long-term (10-year) objectives.

Research Techniques Made Simple: Lipidomic Analysis in Skin Research.

Although lipids are crucial molecules for cell structure, metabolism, and signaling in most organs, they have additional specific functions in the skin. Lipids are required for the maintenance and regulation of the epidermal barrier, physical properties of the skin, and defense against microbes. Analysis of the lipidome-the totality of lipids-is of similar complexity to those of proteomics or other omics, with lipid structures ranging from simple, linear, to highly complex structures. In addition, the ordering and chemical modifications of lipids have consequences on their biological function, especially in the skin. Recent advances in analytic capability (usually with mass spectrometry), bioinformatic processing, and integration with other dermatological big data have allowed researchers to increasingly understand the roles of specific lipid species in skin biology. In this paper, we review the techniques used to analyze skin lipidomics and epilipidomics.

An interferon-related signature characterizes the whole blood transcriptome profile of insulin-resistant individuals-the CODAM study.

BACKGROUND: Worldwide, the prevalence of obesity and insulin resistance has grown dramatically. Gene expression profiling in blood represents a powerful means to explore disease pathogenesis, but the potential impact of inter-individual differences in a cell-type profile is not always taken into account. The objective of this project was to investigate the whole blood transcriptome profile of insulin-resistant as compared to insulin-sensitive individuals independent of inter-individual differences in white blood cell profile. RESULTS: We report a 3% higher relative amount of monocytes in the insulin-resistant individuals. Furthermore, independent of their white blood cell profile, insulin-resistant participants had (i) higher expression of interferon-stimulated genes and (ii) lower expression of genes involved in cellular differentiation and remodeling of the actin cytoskeleton. CONCLUSIONS: We present an approach to investigate the whole blood transcriptome of insulin-resistant individuals, independent of their DNA methylation-derived white blood cell profile. An interferon-related signature characterizes the whole blood transcriptome profile of the insulin-resistant individuals, independent of their white blood cell profile. The observed signature indicates increased systemic inflammation possibly due to an innate immune response and whole-body insulin resistance, which can be a cause or a consequence of insulin resistance. Altered gene expression in specific organs may be reflected in whole blood; hence, our results may reflect obesity and/or insulin resistance-related organ dysfunction in the insulin-resistant individuals.

Comparison of metabolic states using genome-scale metabolic models.

Genome-scale metabolic models (GEMs) are comprehensive knowledge bases of cellular metabolism and serve as mathematical tools for studying biological phenotypes and metabolic states or conditions in various organisms and cell types. Given the sheer size and complexity of human metabolism, selecting parameters for existing analysis methods such as metabolic objective functions and model constraints is not straightforward in human GEMs. In particular, comparing several conditions in large GEMs to identify condition- or disease-specific metabolic features is challenging. In this study, we showcase a scalable, model-driven approach for an in-depth investigation and comparison of metabolic states in large GEMs which enables identifying the underlying functional differences. Using a combination of flux space sampling and network analysis, our approach enables extraction and visualisation of metabolically distinct network modules. Importantly, it does not rely on known or assumed objective functions. We apply this novel approach to extract the biochemical differences in adipocytes arising due to unlimited vs blocked uptake of branched-chain amino acids (BCAAs, considered as biomarkers in obesity) using a human adipocyte GEM (iAdipocytes1809). The biological significance of our approach is corroborated by literature reports confirming our identified metabolic processes (TCA cycle and Fatty acid metabolism) to be functionally related to BCAA metabolism. Additionally, our analysis predicts a specific altered uptake and secretion profile indicating a compensation for the unavailability of BCAAs. Taken together, our approach facilitates determining functional differences between any metabolic conditions of interest by offering a versatile platform for analysing and comparing flux spaces of large metabolic networks.

Investigating the Molecular Processes behind the Cell-Specific Toxicity Response to Titanium Dioxide Nanobelts.

Some engineered nanomaterials incite toxicological effects, but the underlying molecular processes are understudied. The varied physicochemical properties cause different initial molecular interactions, complicating toxicological predictions. Gene expression data allow us to study the responses of genes and biological processes. Overrepresentation analysis identifies enriched biological processes using the experimental data but prompts broad results instead of detailed toxicological processes. We demonstrate a targeted filtering approach to compare public gene expression data for low and high exposure on three cell lines to titanium dioxide nanobelts. Our workflow finds cell and concentration-specific changes in affected pathways linked to four Gene Ontology terms (apoptosis, inflammation, DNA damage, and oxidative stress) to select pathways with a clear toxicity focus. We saw more differentially expressed genes at higher exposure, but our analysis identifies clear differences between the cell lines in affected processes. Colorectal adenocarcinoma cells showed resilience to both concentrations. Small airway epithelial cells displayed a cytotoxic response to the high concentration, but not as strongly as monocytic-like cells. The pathway-gene networks highlighted the gene overlap between altered toxicity-related pathways. The automated workflow is flexible and can focus on other biological processes by selecting other GO terms.

Integrative multiomics analysis of human atherosclerosis reveals a serum response factor-driven network associated with intraplaque hemorrhage.

BACKGROUND: While single-omics analyses on human atherosclerotic plaque have been very useful to map stage- or disease-related differences in expression, they only partly capture the array of changes in this tissue and suffer from scale-intrinsic limitations. In order to better identify processes associated with intraplaque hemorrhage and plaque instability, we therefore combined multiple omics into an integrated model. METHODS: In this study, we compared protein and gene makeup of low- versus high-risk atherosclerotic lesion segments from carotid endarterectomy patients, as judged from the absence or presence of intraplaque hemorrhage, respectively. Transcriptomic, proteomic, and peptidomic data of this plaque cohort were aggregated and analyzed by DIABLO, an integrative multivariate classification and feature selection method. RESULTS: We identified a protein-gene associated multiomics model able to segregate stable, nonhemorrhaged from vulnerable, hemorrhaged lesions at high predictive performance (AUC >0.95). The dominant component of this model correlated with αSMA- PDGFRα+ fibroblast-like cell content (p = 2.4E-05) and Arg1+ macrophage content (p = 2.2E-04) and was driven by serum response factor (SRF), possibly in a megakaryoblastic leukemia-1/2 (MKL1/2) dependent manner. Gene set overrepresentation analysis on the selected key features of this model pointed to a clear cardiovascular disease signature, with overrepresentation of extracellular matrix synthesis and organization, focal adhesion, and cholesterol metabolism terms, suggestive of the model's relevance for the plaque vulnerability. Finally, we were able to corroborate the predictive power of the selected features in several independent mRNA and proteomic plaque cohorts. CONCLUSIONS: In conclusion, our integrative omics study has identified an intraplaque hemorrhage-associated cardiovascular signature that provides excellent stratification of low- from high-risk carotid artery plaques in several independent cohorts. Further study revealed suppression of an SRF-regulated disease network, controlling lesion stability, in vulnerable plaque, which can serve as a scaffold for the design of targeted intervention in plaque destabilization.

Assessing the Contribution of Relative Macrophage Frequencies to Subcutaneous Adipose Tissue.

Background: Macrophages play an important role in regulating adipose tissue function, while their frequencies in adipose tissue vary between individuals. Adipose tissue infiltration by high frequencies of macrophages has been linked to changes in adipokine levels and low-grade inflammation, frequently associated with the progression of obesity. The objective of this project was to assess the contribution of relative macrophage frequencies to the overall subcutaneous adipose tissue gene expression using publicly available datasets. Methods: Seven publicly available microarray gene expression datasets from human subcutaneous adipose tissue biopsies (n = 519) were used together with TissueDecoder to determine the adipose tissue cell-type composition of each sample. We divided the subjects in four groups based on their relative macrophage frequencies. Differential gene expression analysis between the high and low relative macrophage frequencies groups was performed, adjusting for sex and study. Finally, biological processes were identified using pathway enrichment and network analysis. Results: We observed lower frequencies of adipocytes and higher frequencies of adipose stem cells in individuals characterized by high macrophage frequencies. We additionally studied whether, within subcutaneous adipose tissue, interindividual differences in the relative frequencies of macrophages were reflected in transcriptional differences in metabolic and inflammatory pathways. Adipose tissue of individuals with high macrophage frequencies had a higher expression of genes involved in complement activation, chemotaxis, focal adhesion, and oxidative stress. Similarly, we observed a lower expression of genes involved in lipid metabolism, fatty acid synthesis, and oxidation and mitochondrial respiration. Conclusion: We present an approach that combines publicly available subcutaneous adipose tissue gene expression datasets with a deconvolution algorithm to calculate subcutaneous adipose tissue cell-type composition. The results showed the expected increased inflammation gene expression profile accompanied by decreased gene expression in pathways related to lipid metabolism and mitochondrial respiration in subcutaneous adipose tissue in individuals characterized by high macrophage frequencies. This approach demonstrates the hidden strength of reusing publicly available data to gain cell-type-specific insights into adipose tissue function.

A resource to explore the discovery of rare diseases and their causative genes.

Here, we describe a dataset with information about monogenic, rare diseases with a known genetic background, supplemented with manually extracted provenance for the disease itself and the discovery of the underlying genetic cause. We assembled a collection of 4166 rare monogenic diseases and linked them to 3163 causative genes, annotated with OMIM and Ensembl identifiers and HGNC symbols. The PubMed identifiers of the scientific publications, which for the first time described the rare diseases, and the publications, which found the genes causing the diseases were added using information from OMIM, PubMed, Wikipedia, whonamedit.com, and Google Scholar. The data are available under CC0 license as spreadsheet and as RDF in a semantic model modified from DisGeNET, and was added to Wikidata. This dataset relies on publicly available data and publications with a PubMed identifier, but by our effort to make the data interoperable and linked, we can now analyse this data. Our analysis revealed the timeline of rare disease and causative gene discovery and links them to developments in methods.